Jill 


I 




ROLE OF CONDUCTIVITY IN 
HYDROMAGNETIC STABILITY 
OF PARALLEL FLOWS 


by Philip R. Nachtsheim 

Lewis Research Center 
Cleveland, Ohio 


and 

Eli Reshotko 

Case Institute of T echnology 
Cleveland, Ohio 



NATIONAL AERONAUTICS AND SPACE ADMINISTRATION - WASHINGTON, D. C. - DECEMBER 1965 


1 


TECH LIBRARY KAFB, NM 



TECH LIBRARY KAFB, NM 


130132 


NASA TN D-3144 


ROLE OF CONDUCTIVITY IN HYDROMAGNETIC 

STABILITY OF PARALLEL FLOWS 

By Philip R. Nachtsheim 

Lewis Research Center 
Cleveland, Ohio 

and 

Eli Reshotko 

Case Institute of Technology 
Cleveland, Ohio 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


For sale by the Clearinghouse for Federal Scientific and Technical Information 
Springfield, Virginia 22151 — Price $1.00 



ROLE OF CONDUCTIVITY IN HYDROMAGNETIC STABILITY OF PARALLEL FLOWS 


by Philip R. Nachtsheim and Eli Reshotko* 
Lewis Research Center 


SUMMARY 

The role of the electrical conductivity in the stability of magnetohydro- 
dynamic channel flow with parallel magnetic field is examined through exact 
numerical solution of the pertinent sixth-order system of disturbance equations 
throughout an extended range of magnetic Reynolds numbers. The results ob- 
tained indicate that the conductivity of the fluid acts as a stabilizing agent , 
as long as it is small, and as a destabilizing agent, if it is large. It is 
concluded that the conductivity reverses its role as a stabilizing agent at a 
magnetic Reynolds number of order unity. Examination of the contributions to a 
disturbance energy balance equation shows that the conductivity acts as a de- 
stabilizing agent by setting up a time-independent Maxwell stress that has the 
same sign as the vorticity of the basic flow. 


INTRODUCTION 

The hydromagnetic stability of laminar flows of electrically conducting 
fluids has been analyzed by many authors. In particular, it has been shown 
that, when a uniform magnetic field is imposed in the direction of the laminar 
flow, the flow is always more stable than in the absence of a magnetic field. 
This alinement of the magnetic field is especially significant and will be the 
configuration examined herein, since the mean motion of the fluid is not af- 
fected by the Imposed magnetic field; therefore, the net effect of the magnetic 
field on the stability of a given velocity distribution can be investigated. 

Previous treatments of this problem have not made clear the role of the 
conductivity of the fluid in the stability phenomenon. Previous investigators 
have considered two extreme limiting cases of the hydromagnetic stability prob- 
lem; namely, the stability of a fluid of very low and very high electrical con- 
ductivity. 

Stuart (ref. l) treated the case of a fluid at a very low conductivity in 
a very strong magnetic field. His results indicated that the magnetic field 
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stabilized the motion, but his treatment of the problem involved several 
simplifying assumptions, primarily the assumed smallness of the magnetic 
Reynolds number. These assumptions precluded any investigation of the effect 
of the conductivity alone on the stability phenomenon. Velikhov (ref. 2), on 
the other hand, treated the case of a fluid of very high electrical conduc- 
tivity. His results showed that the minimum critical Reynolds numbers were 
markedly greater than at the low values of the conductivity examined by Stuart. 
This would lead one to believe that the stabilizing effect of the magnetic 
field is enhanced as the conductivity is increased. However, Velikhov found 
this trend reversed at the large values of the conductivity investigated; that 
is, an increase at these high values of the conductivity was destabilizing, 
and the minimum critical Reynolds number for infinite magnetic Reynolds number 
is below that obtained for large but finite magnetic Reynolds number. However, 
with Velikohov's simplifying assumptions, primarily the assumed largeness of 
the magnetic Reynolds number, again no indication could be obtained of when 
the conductivity reversed its role as a stabilizing agent. 

In the present report, the role of the conductivity is clarified by ex- 
amining the stability problem between the limiting cases of very low and very 
high conductivity. This is accomplished by solving the complete disturbance 
equations numerically without invoking the simplifying assumptions made in 
earlier treatments of this problem. 


SYMBOLS 

Vector components with a numerical index refer to dimensional quantities; 
vector components with a literal index refer to nondimensional quantities. 

A Alfven number, yB 2 /pu^ 0 

B reference magnetic induction 

b-pbg disturbance magnetic induction 

c c r + ic^ 

Cj_ time amplification factor of disturbance wave 

c r phase speed 

E electric field 

£ disturbance energy 

J electric current density 

k wave number 

L reference length 
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N magnetic interaction parameter , cB^L/pl^ 

P pressure 

Pr m magnetic Brandtl number, o\i Q v 

p disturbance pressure 

Re Reynolds number, U^L/v 

Re m magnetic Reynolds number, 

o 

s disturbance vorticity amplitude, cp Tt - a cp 

t time 

U velocity 

IJjjjL reference velocity 

u-^,Ug disturbance velocity components 

Xi,x 2 Cartesian coordinates, dimensional 

x,y Cartesian coordinates, nondime nsional 

a dimensionless wave number, Lk]_ 

£ disturbance vorticity 

0 disturbance electric current density amplitude, 

p 0 permeability (vacuum value) 

v fluid kinematic viscosity 

IT P + (B • B)/2n 0 

it disturbance total pressure amplitude 

p fluid density 

a electrical conductivity 

$ disturbance stream function 

cp disturbance stream function amplitude 

1 disturbance vector potential 

\| ; disturbance vector potential amplitude 
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Subscripts : 


i imaginary part 

r real part 

c differentiation with respect to c 

cr critical 


FORMULATION OF PROBLEM 


The general equations of magnetohydrodynamics are the equations of elec- 
trodynamics for moving media and the Navier-Stokes equations modified to in- 
clude the electromagnetic body force. For an incompressible fluid with a 
scalar conductivity the governing equations are (ref, 3) 



+ U 


vuj =-VP + JxB + pv 

CD 

v x E " " sf 

(2) 

v x H = (i 0 J 

(3) 

J = o(E + IJ X B) 

(4) 

o 

II 

Tim 

D> 

(5) 

V • B = 0 

(6) 


The kinematic viscosity coefficient v is taken to be a constant, and the ex- 
cess charge density and displacement currents are neglected in the preceding 
equations in accordance with the usual magnetohydrodynamic approximations. 


The steady flow under consideration is the flow with parabolic velocity 
profile between perfectly conducting parallel planes in the presence of a con- 
stant Imposed magnetic field parallel to the fluid velocity. 

The solution of the steady equations that satisfies the imposed conditions 
is 


U : [u 1 (x 2 ), 0, o] 

5 : [B-l, 0, 0] 

where is the constant imposed magnetic field and 
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( 7 ) 



where 2L is the spacing between the planes and Ufo is the velocity at the 
centerline (x 2 = 0). It is now proposed to consider the equations governing 
small two-dimensional disturbances u : (uq,u 2 ,0), Id ; (bq,b 2 ,0), and p, which 
are superimposed on the steady-state solution. There is no loss in generality 
in assuming that the disturbances are two-dimensional in nature , since , for the 
problem under consideration, the motion is always more stable for three- 
dimensional than for two-dimensional disturbances as pointed out by Stuart 
(ref. l). This is the analogue of Squire's theorem (ref. 4) established in the 
ordinary incompressible parallel flow stability theory. Periodic sinusoidal 
disturbances are considered in which u, b, and p are all functions of X 2 
multiplied by exp[jLkq(xq - Cqt )j , where cq is the complex phase velocity and 
is the wave number in the x^-direction. Since kq is always positive, 
the disturbances are amplified or damped according to whether the imaginary 
part of Cq is positive or negative. If the Imaginary part of cq is zero, 
the disturbance is neutrally stable. The steps leading to the disturbance 
equations are carried out in the appendix. This procedure consists of super- 
imposing u, b, and p onto the steady-state solution, substituting them into 
the general equations (l) to (6), and linearizing the equations with respect to 
the small disturbances It, b, and p. Elimination of the variables uq, bq, 
and p leads to two simultaneous ordinary differential equations for b 2 and 
u 2 . If U m cp(x 2 ) a ^d Bq\|/(x 2 ) denote quantities proportional to the amplitude 
functions of u 2 and b2, respectively, the following nondimensional simulta- 
neous equations are found to govern cp and ^ 

cp”” - 2aV + a 4 cp = ia£e[(u x - c)(cp" - a 2 cp) - U^cp - A 2 (f" - a 2 i|0] (8) 

V - a 2 t = iaRe m j(U x - c)i|r - cpj (9) 

where x = Xq/L, y = x 2 /L, c = c-^/u^ a = Lkq, Re = l^L/v is the Reynolds 

number, Re m = a^U^L is the magnetic Reynolds number, and A 2 = B q/pU^p 0 is 

the Alfven number squared. The primes denote differentiation with respect to 
y. It is to be noted here that, in the absence of a magnetic field (A* = 0), 
equation (8) reduces to the ordinary Orr-Sommerfeld equation. 

Two new parameters enter into the consideration of the stability of the 
flow of an electrically conducting fluid; namely, the magnetic Reynolds number, 
which is proportional to the conductivity of the fluid, and the Alfven number, 

which is proportional to the strength of the imposed magnetic field. 

Equations (8) and (9) and the boundary conditions (appropriate for per- 
fectly conducting planes which bound the flow) 

y = ±l : cp=cp ! = \|r = 0 (10) 

constitute an eigenvalue problem, which, for a, Re, Re m , and A fixed, con- 
sists of determining c in order to satisfy the boundary conditions. 

Earlier treatments of the eigenvalue problem have considered the limiting 
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cases Re m « 1 (ref. l) or Re m » 1 (ref. 2). In both limiting cases, the 
governing equations can be simplified somewhat. The resulting simplified equa- 
tions are more amenable to solution than the coupled equations (8) and (9). 
However, if it is proposed to examine the behavior of the solutions of equa- 
tions (8) and (9) between the limiting values of Re m « 1 and Re m » 1, the 
complete set of equations must be considered. 


Simplified Disturbance Equations 

The assumption underlying the simplified treatment of the eigenvalue prob- 
lem in reference 1 will now be examined. For many electrically conducting 
fluids or slightly ionized gases used in laboratory experiments, Re m is 
usually very small. However, the assumption Re m « 1 is by itself not suf- 
ficient to derive Stuart's form of the disturbance equations. Stuart's simpli- 
fied disturbance equation can be obtained from equations (8) and (9) by fol- 
lowing the method of Tatsumi (ref. 3); namely, elimination of cp, making the 
transformation 6 = ty" - a/ty, and performing the limiting process Re m -* 0 
but A 2 Re m finite. These steps result in the following equation for 0: 

0"" - 2a 2 0” + a 4 0 = iccRe [(U x - c)(0" - a 2 0) - U x 0 + iaA 2 Re m 0] (ll) 

The boundary conditions for equation (ll) follow from equations (8) and (9) 
and equation (10) if the limiting form of these equations is considered. The 
boundary conditions are 


y = ±1 : Q - Q x =0 (12) 

The transformation and subsequent limiting process clearly show the nature of 
the assumptions underlying this simpler eigenvalue problem. In terms of physi- 
cal properties, the assumptions are the flow of a fluid of low conductivity and, 
in addition, under a very strong magnetic field. It is to be noted here that 
the quantity A 2 Re m = aB^L/pl^ = N is referred to as the magnetic interaction 
parameter. When this parameter is negligibly small, equation (ll) is identical 
to the Qrr - Sommer f eld equation. Other properties of this extreme limiting case 
are that the two parameters A 2 and Re m have been collapsed into a single 
parameter, N = A 2 Re m , and that the order of the coupled differential system 
(eqs. (8) and (9)) has been reduced from six to four. Also the six boundary 
conditions (eq. (lO)) have been used in conjunction with equation (9) to formu- 
late the four boundary conditions (eq. (12)) for the case of the infinitely 
conducting wall. 


Numerical Solution of Eigenvalue Problem 

The numerical solution of the eigenvalue problems formulated will now be 
considered. The general method used is quite analogous to the method estab- 
lished in reference 5 for solving the Orr - Sommer f eld equation; therefore, it 
will be sufficient to give only a brief account of the steps required to gen- 
eralize the established method in order to handle this problem. Furthermore, 
since equation (ll) is a limiting form of equations (8) and (9), the present 
account will be concerned with the more general system (eqs. (8) and (9)) 
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together with the boundary conditions of equation (lO). 


Before describing the numerical method, a simplification of the problem, 
which halves the range of integration, should be pointed out. This simplifica- 
tion is based on the observation that the disturbances in equations (8) and (9) 
can be separated into even and odd functions, since the variable coefficients 
in these equations, namely U x and U x , are even functions of y. In ordinary 
hydrodynamic stability theory, the even solution, which has the simpler flow 
pattern, usually gives a lower minimum critical Reynolds number and will be the 
solution examined herein. This simplification enables one to consider only 
even solutions in half of the channel through introduction of appropriate sym- 
metry conditions at y = 0. The new boundary conditions for equations (8) and 
(9) are 


y = 0 : cp 1 = cp 111 = ^ 1 =0 


(13) 


y = 1 : cp = cp 1 = \|r = 0 (14) 

2 

The approach to the eigenvalue problem for fixed a, Re, A , and Re m 
used herein is to find values of c = c r + ic^ (eigenvalues) for which equa- 
tions (8) and (9) have solutions (eigenfunctions) that satisfy the boundary 
conditions equations (13) and (14). 


In order to find an eigenvalue, the following iterative procedure is car- 
ried out. A trial solution of equations (8) and (9) is obtained by numerical 
integration. In order to integrate numerically, additional boundary values are 
specified at both boundaries y = 0 and y = 1, and a value of c is speci- 
fied. Then both integrals are stepped in toward the middle y = l/2. The addi- 
tional boundary values and the eigenvalue c have to be adjusted in order to 
match the solutions in the middle. After a trial integration, subsequent 
boundary values and eigenvalues are automatically calculated by the Newton- 
Raphson technique of finding successive approximations. 


The various elements that enter into the preceding iterative procedure will 
now be put forward. Instead of solving equations (8) and (9) as they stand, a 
system of second-order differential equations is formulated by introducing the 
disturbance vorticity amplitude function s = cp" " a,^cp. In terms of the vorti- 
city amplitude, the system (eqs. (8) and (9)) goes over to 


cp" = a 2 cp + s 


\| r" = a 2 f + iccRe m ^(U x - c)i - cp] 
s" = a 2 s + iaRe [(U x - c)s - U^cp - A 2 (\|r" - a 2 ^)] 


(15) 

( 16 ) 
(17) 


In order to obtain a system of equations in which derivatives of the dependent 
variables do not appear on the right side, \|r" is eliminated from equation (17) 
by means of equation (16) giving 

s" = a 2 s + iaRe[(U x - c)(s - iaNt) - (U" - iaW)cp] (18) 
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where 


N = A 2 Re m 

Equations (15), (16), and (18) are solved subject to the boundary conditions 
(eqs . (13) and (14) ) that employ the definition of the vorticity amplitude 
function. 

For a forward solution (y increasing) starting at y = 0, the initial 
values are specified according to the following table 

y = 0 : 

For a backward solution (y decreasing) starting at y = 1, the initial 
values are 


y = 1 

The additional boundary values that have to be specified in order to integrate 
numerically are those enumerated in the second row of each table. The condi- 
tion cp(0) = 1 is a normalizing condition and fixes the size of the whole so- 
lution. With the preceding boundary values specified and a value of c speci- 
fied, trial solutions of equations (15), (16), and (18) are obtained. Next, 
the process of matching at a common point is carried out. The matching point 
was taken to be y = l/2 to equalize the numerical errors that grow in propor- 
tion to the number of integration steps taken. At the matching point, the so- 
lution must be continuous. The requirements of continuity lead to six condi- 
tions that have to be satisfied at y m l/2. The six quantities S(0), ty(0), 
s(l), s'(l), \|r r (l), and c have to be adjusted in order to satisfy these con- 
ditions. Let the difference between the forward and backward solutions eval- 
uated at y = l/2 be designated by gj_(i = 1 ... 6) and set up a correspon- 
dence between the differences and the various functions according to the fol- 
lowing table: 


S6 

Now, the gj_ depend on the adjustable parameters (the additional boundary con- 
ditions and c). The Newton-Raphson method is used to fulfill the condition of 
continuity at the matching point (i.e., all g^ = 0). 

For ease in writing, the various adjustable parameters are designated ac- 
cording to the following scheme: 
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s(0) 

♦(0) 

s(l) 

s'(l) 

t'(l) 

c 

Z 1 

Z 2 

Z 3 

Z 4 

Z 5 

Z 6 


If the chosen values z^ produce a solution that gives all g. as approxi- 
mately zero, a better approximation is obtained by starting "with + Azq. 
The quantities Azq are solutions of the equations 



(i = 1 . . .6) 


( 19 ) 


The partial derivatives are obtained by solving additional initial- 

value problems. All the partial derivatives with respect to the additional 
boundary conditions z^(i = 1 . . .5) can be obtained by solving equations 
(15), (16 ) and (18) with appropriate initial conditions since these equations 
are linear. For example, derivatives with respect to s(0) of the variables in 
the forward integration are obtained by solving equations (15 ) , (16), and (18) 
starting at y = 0 with initial conditions given according to the following 
table : 


cp 

cp' 

s 

s 1 


u 

0 

0 

1 

0 

0 

0 


For the backward integration, the variables are taken to be independent of 
s(0). For derivatives with respect to c, a different system of equations has 
to be solved. This system is obtained by partial differentiation of the terms 
in equations (15), (16), and (18) with respect to c. The quantities 



satisfy the system of equations 

cpc = cc 2 cp c + s c (20) 

= + iaRe mP U x “ c )t c ~ 9 C “ *)] (21) 

s c = a2s c + iaJRe [( u x " c )( s c " iaWtc) " (U" - i<xN)cp c - (s - iaN\|r)] (22) 

The appropriate initial conditions for equations (20) to (22) are homogeneous 
initial conditions for both the forward and backward integrations. 

After obtaining a trial solution and the various partial derivatives, 
equation (19) can be solved for the corrections. Successive applications of 
this procedure should converge to an eigenvalue c. 

p 

The procedure just outlined for finding an eigenvalue c for a, Re, A , 
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11 000 , 


and Re m fixed -was programed for solu- 
tion on the IBM 7094 computer located 
at the Lewis Research Center. 

In reference 5, the question of 
the accuracy and rate of convergence of 
the method is discussed with regard to 
solutions of the Orr-Sommerf eld equa- 
tion. Also, several simplifications 
are pointed out there that reduce the 
amount of labor that is required in the 
application of the method. 

Carrying out the foregoing proce- 
dure for a system of second-order 
equations enables one to obtain exact 
numerical solutions of equations (8) 
and (9), and equation (ll). 


SOLUTIONS 

Solutions for Re m « 1 

A comparison will now be made of 
the exact numerical solutions of equations (8) and (9), the sixth-order sys- 
tem; the exact numerical solutions of equation (ll), the fourth-order system; 
and Stuart ! s asymptotic solution of the fourth-order system. This comparison 
is made under the conditions that should satisfy the assumptions made by Stuart 
in going from the sixth-order system to the fourth-order system. These condi- 
tions would correspond to the flow of mercury with a magnetic Prandtl number 
Rr m of 1.5xl0~ 7 . For a Reynolds number of 10 4 , the value of the magnetic 
Reynolds number would be about 1.5x10“ 3. 

Figure 1 shows the minimum critical Reynolds number as a function of the 
interaction parameter N corresponding to the solutions of the three eigenvalue 
problems. The significance of the minimum critical Reynolds number is that all 
disturbances will be damped in flows with Reynolds numbers below the minimum 
critical value. The exact numerical solutions of the sixth- and fourth-order 
systems yield identical results. Thus, it can be inferred that Stuart*s trunca- 
tion from the sixth-order to the fourth-order system is justified when 
Re m 1 • 

Also shown in figure 1 are Stuart ! s results for this case. The difference 
in the values of Re cr at the same value of the interaction parameter is due 
only to the difference in the methods used to solve equation (ll). The results 
obtained by Stuart differ from those obtained by the present method by as much 
as 30 percent. Furthermore, the value of Re cr obtained by the present method 
for N = 0, in which case equation (ll) reduces to the Orr-Sommerf eld equation, 
agrees with the value of Re cr , which Thomas (ref. 6) obtained by solving the 
Orr-Sommerf eld equation. 



Figure 1. - Minimum critical Reynolds number. 
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Figure 2. - Neutral stability diagrams for very large magnetic Reynolds numbers (from 
ref. 2). Alfven number, 0.08. 


Solutions for Re m Not Small 

The solution of the complete eigenvalue problem corresponding to the dif- 
ferential equations (8) and (9) (i.e., the sixth-order system) and boundary 
conditions will now be discussed for the case where Re m is not restricted to 
be small. The case Re m » 1 was solved by Velikhov by obtaining asymptotic 
solutions of equations (8) and (9) at large values of aRe and ocRe m - The re- 
sults obtained by Velikhov for an Alfven number, A = 0.08, are shown in fig- 
ure 2. As may be seen, the minimum critical Reynolds numbers are considerably 
higher than those for Re m « 1 shown in figure 1. Hence, it could be con- 
cluded that the stabilizing effect is enhanced as Re m is increased. However, 
as can be seen from figure 2, the higher Re m yield the lower Re cr . As 
pointed out by Tatsumi (ref. 3), it seems probable that there exists some value 
of Re m above which the role of conductivity as a stabilizing agent is re- 
versed. From the various factors considered subsequently and the results of 
the next section on the balance of disturbance energy, it appears that the con- 
ductivity reverses its role as a stabilizing agent at values of Re m of order 
unity. 

The question of reversal of the stabilizing effect of conductivity is ex- 
amined by fixing Re and a and determining the variation of c^ for in- 
creasing Re m in numerical solutions of equations (8) and (9). A destabiliz- 
ing effect of increasing Re m will be evidenced by Cj_ increasing. This 
would be equivalent to Velikhov* s calculated decrease of Re cr as Re m is in- 
creased. Velikhov, however, made his calculations after passing to the limit 
of large Re m and hence could not give an indication of when the conductivity 
reversed its role as a stabilizing agent. In the present study, on the other 
hand, the entire range of Re m from Re m « 1 to Re m » 1 has been investi- 
gated by varying the magnetic Hrandtl number. 

The results of the present calculations are shown in figure 3. These cal- 
culations were carried out at a fixed point in the a, Re plane; that is, 
a = 1, Re = 9000. Since Re was held constant, each value of Re m shown in 
figure 3 corresponds to a different magnetic Prandtl number, Pr m 
(Pr m = Re m /Re = pcrv). Scales of both quantities Re m and Pr m are shown 
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Figure 3. - Effect of increasing magnetic Reynolds number. Nondimensional wave number, 1; Reynolds number, 9000. 


as abscissae in figure 3, where c^, the time amplification factor is shown as 
the ordinate. Curves are shown both for constant values of the magnetic inter- 
action parameter N, and for the Alfven number A, which is directly propor- 
tional to the magnetic field strength. 

The results presented in figure 3 may be interpreted in the following 
manner: Imagine, for example, a flow in a given channel where the fluid vis- 

cosity and flow velocity are fixed and the flow is subjected to disturbances of 
a given wavelength. This corresponds to a fixed point in the a. Re plane. 

In this example, the conductivity of the fluid and the strength of the Imposed 
magnetic field can be altered at will. The point a = 1, Re = 9000 is an 
unstable point in the absence of a magnetic field, as shown by the curve A = 0, 
N = 0 in figure 3. 

To be noted immediately is that, for an electrically conducting fluid, the 
presence of any magnetic field tends to reduce the amplification rate below 
that obtained in the absence of a magnetic field. This is true regardless of 
whether the magnetic Reynolds number is large or small. For a given magnetic 
Reynolds number, the damping rate relative to the nonmagnetic case increases 
monotonically with an increase in magnetic field strength; for a given magnetic 
field strength (represented by constant A in our example) the damping rate 
increases monotonically with increasing interaction parameter. The magnitudes 
of these variations, however, depend on whether the magnetic Reynolds number is 
small or large. 
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For a constant interaction parameter N, the variation of amplification 
rate with magnetic Reynolds number is negligible until Re m increases 

past order unity, whence it increases sharply relative to the value for 
Re m « 1 and gradually approaches the value in the nonmagnetic case as Re lU 
becomes very large. The present results for Re^ « 1 are consistent with 
those obtained from the fourth-order system (eq. (ll)). It is seen that cp 
depends only on the single parameter N when Re m « 1, but varies also with 
Re m when Re m exceeds order unity. As Re m increases at constant N, the 
growth rates increase, but are always less than in the nonmagnetic case. 

For constant A, which in our example is interpreted as constant magnetic 
field, increasing the conductivity always gives damping. However, there is 
substantial increase in the damping rate with the magnetic Reynolds number when 
Re m is small; for large magnetic Reynolds number, an increase in Rem tends 
to increase the damping only slightly. 

A more detailed physical description of the nature of the stability prob- 
lem as Re m increases past unity may be obtained by examining the terms in a 
disturbance energy balance. This examination is made subsequently. 

The general situation here in which two transport properties enter, 
namely, the electrical conductivity and the viscosity, is quite analogous to 
the situation in which the viscosity coefficient alone enters. In ordinary 
hydrodynamic stability, it is known that viscous mechanisms are stabilizing 
(the dissipation) as long as the viscosity is large enough (Re small enough). 
It is only when the viscosity becomes small enough (Re large enough) that 
viscous mechanisms come into play which destabilize the resulting motion by 
providing a means of transferring energy to the disturbance motion through the 
Reynolds stress. In magnetohydrodynamic stability, the role played by the re- 
sistivity (the reciprocal of the conductivity) is analogous to the role played 
by the viscosity in ordinary hydrodynamic stability; that is, at small Re m 
(high resistivity) the conductivity o is primarily stabilizing because of the 
predominance of the joule heating. At large Re m , the conductivity tends to be 
destabilizing under certain conditions because there is a mechanism that can 
transfer energy to the disturbance motion. 


DISTURBANCE ENERGY BALANCE 

Formulation of a disturbance energy balance and subsequent examination of 
the various contributions to the rate of change of disturbance energy gives an 
insight into the mechanics of the disturbance motion and enables one to identify 
the various processes that are operating. 

The starting point for the derivation of the disturbance energy balance 
equation are equations (A9) to (A12). Multiplying the terms in equations (A9) 
by u xj (A10) by Uy, (All) by A 2 b x , and (A12) by A 2 by, adding, and perform- 
ing some algebraic manipulations yield (for B x = l) : 


13 


[ 


& e - - K'V - A=V d £-i;l z - 
- s (u ^’ t) " 4 ( V> + I? [s (u y ?) 
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where 

u§ + u? 

= ^ 3 "y , 2 D x _ p y 
^ _ 2 + A 2 


and 

Dg _ Be Be 
M » 3t + Ox £ 


The quantity G is the dimensionless kinetic energy of the fluid plus the 
energy of the magnetic field, all referred to the kinetic energy of the mean 
flow. The quantity £ = duy/dx - chi x /c)y is the dimensionless vorticity of the 
disturbance flow, and j = oby/^x - db x /dy is the dimensionless curl of the 
disturbance magnetic field, which by Amperes law is equal to the disturbance 
electrical current density. 


Equation (23) gives the time rate of increase of the disturbance energy 
(per unit volume) of a fluid element that moves with the basic flow. The terms 
in this equation are to be integrated over a cell that extends across the chan- 
nel in the y-direction and along the channel for a distance of wavelength. 

All terms in the braces in equation (23) vanish either because a disturbance 
quantity vanishes at the boundaries of the flow for the integration across the 
channel or because of periodicity in the direction along the channel. The 
equation giving the time rate of increase of energy of the disturbance motion 
is therefore : 




dx dy 


f dU- 

" u x u v 

x y dy 


x 


o AU* 

+ A b x by — 


-L # 

Re ^ 



(24) 


Of course, it is the real part of each disturbance amplitude that is required 
in equation (24). 


It is appropriate at this time to identify the various factors on the 
right side of equation (24), which contribute 'to the rate of change of the dis- 
turbance energy. As in ordinary hydrodynamic stability theory, there appear 
the Reynolds stress term and the viscous dissipation term, the first and third 
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terms, respectively, on the right side of equation (24). As is well known, 
the dissipation term is always stabilizing, and the Reynolds stress term can "be 
either stabilizing or destabilizing. The new terms that appear in the stability 
of a conducting fluid are the second and fourth terms. The fourth term 
(-A 2 j 2 /Re m ) can be identified as joule heating and appears in the energy balance 
equation in such a way that this term mil always be stabilizing. The quantity 
b x by in the second term can be identified as a component of the Maxwell stress 
tensor of electromagnetic theory. The question of whether this new term is 
stabilizing or destabilizing can be answered by solving the eigenvalue problem 
and computing the eigenfunctions. 


The disturbance quantities on the right side of equation (24) can be ex- 
pressed in terms of the complex amplitude functions cp and ijr by means of the 
relations 


u x = (y)exp^ia(x - ct)j| 

Uy - RZ'|-iacp(y)expj^ia(x - ct)j| 
b x = Rl|^' (y)expjOa(x - ct)]j 
by = Rl |-iai!r(y)exp jji.a(x - ct)j| 


(25) 

(26) 

(27) 

( 28 ) 


where HI denotes the real part of the complex quantity. The integration with 
respect to x can be performed in equation (24) after substitution of equa- 
tions (25) to (28). For neutral disturbances (c^ = 0) 



9 dU. 

+ A^(\lr , i)/. - \|/ ! ) n 

^ Y r Y i hr h' dy 


x 


aRe 


( 9 . 


r 


a 2 cp r )' 


(TV - 


a 2 cpp 


Ar 


aRe 




iH t T + (q - a 2 *.)' 


(29) 


The range of integration is taken from y = 0 to 1 since all the integrands 
are even functions. For other than neutral disturbances } it is only necessary 
to multiply the right side of equation (29) by factor exp(2ac-jt) in order to 
obtain the time rate of change of disturbance energy. 


In view of the conclusion drawn in the previous section concerning the role 
played by the conductivity as a stabilizing agent , it is of interest to com- 
pare the various terms that contribute to the disturbance energy balance when 
Re m is small and when it is large. Figure 4 shows the distribution across 
half the channel of the electromagnetic terms in the energy balance equation for 
the two cases Re m = 0.00135 and 88.5 (N = 0.0216 for both). The area under 
these curves represents the net contribution to the rate of increase of the 
disturbance energy. For Re m = 0. 00135 , the Maxwell stress is negligible and 
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Normal distance from center of channel, y 

Figure 4. - Maxwell stress and Joule dissipation. 
Nondimensional wave number, 1; Reynolds 
number, 9000; magnetic interaction parameter, 
0.0216. 



Figure 5. - Reynolds stress and viscous dissipation. 
Nondimensional wave number, 1; Reynolds num- 
ber, 9000; magnetic interaction parameter, 0.0216. 


is not plotted, but the joule dissipation contributes substantially to the 
stabilization of the motion (the area under the curve being negative in this 
case). For Re m = 88.5, the Maxwell stress becomes significant and is de- 
stabilizing. Of course, the joule heating is again stabilizing, but the con- 
tribution is much less than in the previous case.. 

Figure 5 shows the distribution across half the channel of the hydrodyna- 
mic terms in the energy balance equation for the two cases. The distribution 
of the viscous dissipation energy is very closely the same for the two Re m 
values, but, for the Reynolds stress energy distributions, there is less area 
under the curve for Re m = 88.5 than for Re m = 0.00135. However, this de- 
crease in production of disturbance energy by the Reynolds stress is offset by 
the destabilizing effects of the decrease in the joule dissipation and the in- 
crease of disturbance energy by the Maxwell stress (fig. 4). The net effect of 
increasing the conductivity is then a consequence of the balance between the 
Reynolds stress mechanism on the one hand and the joule dissipation and Maxwell 
stress mechanism on the other hand. This net effect can be determined by esti- 
mating the areas under the curves in figures 4 and 5. This comparison leads to 
the result that, at a fixed point in the a, Re plane, the net effect of in- 
creasing the conductivity N (the magnetic interaction parameter) constant is 
destabilizing in agreement with the results presented in figure 3. As the Re m 
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is raised.^ the "behavior of the various terms in the energy balance equation 
under the conditions mentioned previously can be summarized as follows: 

Reynolds stress: stabilizing in that there is less production of disturbance 

energy; viscous dissipation: negligible change; Maxwell stress: destabilizing; 

production of disturbance energy occurs only at high values of Re m ; joule dis- 
sipation: destabilizing; the amount of energy that can be dissipated is re- 
duced as the conductivity is raised. If N is not held constant but is allowed 
to increase as Re m increases; the flow is again stabilized as observed by 
Velikov (ref. 2). 


CONCLUSIONS 

The stability of plane magnetohydrodynamic channel flow with parallel mag- 
netic field has been reexamined through exact numerical integration of the 
pertinent sixth-order system of disturbance equations; subject to appropriate 
boundary conditions. The results confirm that Stuart’s reduction of the prob- 
lem to a fourth-order disturbance equation is valid for magnetic Reynolds num- 
ber small compared with 1. However; Stuart’s asymptotic values are about 
30 percent below the present numerical results. For magnetic Reynolds numbers 
of order 1 or greater; there are significant changes in the stability charac- 
teristics. This is borne out by a calculation of the various viscous and mag- 
netic' contributions to the rate of change of disturbance energy. 

It is shown that the resistivity enters this problem in two ways: (l) it 

sets up a time -independent Maxwell stress that augments the disturbance energy 
when it is of the same sign as the vorticity of the basic flow; and (2) through 
joule dissipation the disturbance energy is decreased. For small resistivity 
(Re m » l) the former dominates; leading to a net augmentation of disturbance 
energy; while for large resistivity (Re m « 1) the dissipative effect is 
dominant. 


Lewis Research Center; 

National Aeronautics and Space Administration; 
Cleveland; Ohio; September 22; 1965. 
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APPENDIX - DERIVATION OF DISTURBANCE EQUATIONS 

Although the disturbance equations used herein appear in the literature 
(e.g., ref. 3), they must be displayed in component form in order to formulate 
a disturbance energy balance equation. Since this form appears in the process 
of deriving the disturbance equations, a brief derivation of these equations 
follows . 


The starting point for the derivation is equations (l) to (6), 
tion of J from equations (l) and (3) results in 


Elimina- 


OU - 

P §t + U 


VU 


-vn + — b 

Bo 


VB + pvV^U 


(Al) 


where 


H = P + 


B • B 
2 Do 


Since two-dimensional disturbances are being considered, only a two-dimensional 
magnetic field need_^be considered. A second equation is obtained by the elimi- 
nation of E and J from equations (2) to (4) with the aid of equations (5) 
and (6) and results in 

— >■ 

H + U-VB = B-W + — ^B (A2) 

at crpo v ' 

The stability investigation is limited to mean flows that satisfy equation (7). 

The disturbance equations are obtained by introducing a two-dimensional 
disturbance into equations (5), (6), (Al), and (A2 ) ; that is, let 

tJ -* 1? + u j ( Uq + Uq , u 2 , 0 ) 

S -> S + b : (®q ^qj) b 2 , t 

P -» P + p 


where the lower-case quantities are small perturbations of the basic quantity. 
After this is accomplished and after subtracting out the basic flow, the dis- 
turbance equations are 


P 



U2 


dU-i 

t — — + U 1 

dx 2 




+ — Bq 
Do 


<db-j- 

SXq 


„ 2 

pVV Uq 


(A3) 



+ V 



V - + — B 1 
Sx 2 D 0 ! 


Sbg 


pvV 2 Ug 


(Ai) 


18 



where 


Sb-i 




ciF" + u i ^ = 13 


dU l . „ Su l . 1 „ 2 


3Xn 


2 d.X2 + Bl Sx 1 


V^bn 

a Rn 1 


ab- 

c 

at 


Sb- 


+ u- 


i ax. 


= B 


Su 2 1 o 

■j "n V bp 

-L OX. 0JJ.Q < r: 


Sx 1 + dx 2 ~ ^ 


abq Sb^ 
dxq + dxg - ^ 


it = p + 


E * b 

P-n 


(A5) 

(A6) 

(A7) 

(A8) 


It is convenient to cast the disturbance equations into nondimens ional form by- 
making the following replacement of variables: 


Xq -» xL y xg -> yL, t 




u- 


Uq “* UxHjx; Bq B X B 

1 - u xV u 2 - u y U m^ b l - b x B ' b 2 V 
* - (pUm)* 

where L, and B are fixed reference dimensional quantities. The dis- 

turbance equations in nondimensional form are 


au. 


dlL 


au. 


^ 'V' ^ V 

5r + u ydT + u x5T 


+ a2-d ^ b x 1 

^ + AB x^r + fte V U * 


SUy 

5t + u x 37 


db. 


X 

ab. 


a^t 

37 


+ a 2 b. 


Sb 


y 


+ Jkr V 2 ’ 


'x 5F" Re V ^y 


dU. 


x TT ~~x , 

5t~ x ax y dy 


x + B x ^E. + v 2 b x 
x dx Re m x 


ab q p 

St“ + u x §3T = + B x + Ri^ v b y 


(A9) 

(A10) 

(All) 

(A12) 
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(A13) 


where 


is the Alfven number squared 


is the Reynolds number, and 


Su x 5u y 

5x dy 


0 


db,, db. 


x 


Sx sr 


= o 


A 2 = 




P U m^o 


Re 




Re m - 0'q o l^ ri L 


is the magnetic Reynolds number. 


(A14) 


Besides these basic nondimensional quantities, certain other combinations 
of them are often used; namely, 

^m = °Po v 

the magnetic Prandtl number and 


N = 


ctB 2 L 

P^n 


the magnetic interaction parameter. The following relations apply: 

Re m = * Re 


and 

N = A 2 Re m 

The number of dependent variables can be reduced by introducing a stream 
function for the disturbance velocity and a vector potential for the magnetic 
field. Let 


u x = 


¥ 




U y " " ^ 


(A15) 

(A16) 
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(A17) 


_ S'? 
x " 

CM 

'y ” " Sx 


(A18) 


Substitution of equations (A15) to (A18) into equations (A9) to (A12) leads 


to the following disturbance equations : 

dlly ^ 9 1 

" dy + U x 0 xy 3x + A B x^yx + Ri" ^yxx + $ yyy) (A 19 ) 

®xt + D x®xx - + 5p + A Vn + g? + ®x»y> C 420 ) 

dU x -i 

Y yt + U x^xy = " ay” + B x®xy + rT^ ^yxx + ^yyy) ( A2 l) 

^xt + ^x^xx = ®x®xx + Re m (^xxx + ^xyy) (A22) 


Elimination of jt between equations (A19) and (A20) by cross differentiations 
leads to 

cL^U 

®yyt + ^xxt + ^x^xyy + ^xxx) “ 2 = [®x(^xyy + ^xxx)] 

^-«y 

4- i (0 4-20 4-0 ) (A23) 

Re v yyyy xxyy ^xxxx ; v * 

If the form of the disturbance is taken to he 

0 = cp(y)exp[ia(x - ct)J (A24) 

¥ = ^(y)exp [ia(x - ct)] (A25) 

equation (A23) becomes 

cp ,,n - 2a 2 cp M 4- a^cp = iccRe|(U x - c)(cp" - a 2 cp) - U x cp - A 2 [b^V 1 - a 2 ^(r)]j' 

(A26) 

where the prhnes denote differentiation with respect to y^ and equations (A2l) 
and (A22) become 

V - a 2 t = ^^[(Ux “ c )+ " (A27) 

Equations (A2l) and (A22) lead to a single independent equation since equa- 
tion (A2l ) , for the assumed form of the disturbance , can be obtained from equa- 
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tion (A22) "by differentiation with respect to y. 

Equations (A26) and (A27) are the final form assumed by the disturbance 
equations. Since is a constant, it is appropriate to replace B x by 

unity; that is, the reference quantity B is taken to be Bj_. This final step 
yields equations (8) and (9). 
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